Measurement of magnetic excitations in the two-dimensional antiferromagnetic 
Sr2Cu02Cl2 insulator using resonant x-ray scattering: Evidence for extended interactions 
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We measured the momentum dependence of magnetic excitations in the model spin- 1/2 2D antiferromagnetic 
insulator Sr2Cu02Cl2. We identify a single-spin- wave feature and a multi-magnon continuum, with different 
polarization dependences. The spin-waves display a large (70 meV) dispersion between the zone-boundary 
points (7r,0) and (7r/2,7r/2). Employing an extended t-t'-t"-U one-band Hubbard model, we find significant 
electronic hopping beyond nearest-neighbor Cu ions, indicative of extended magnetic interactions. The spectral 
line shape at (7r,0) indicates sizeable quantum effects in SCOC and probably more generally in the cuprates. 



Magnetism in low-dimensional cuprates remains of con- 
siderable interest, in relation both to the fundamental quest 
to understand strong electron correlation and quantum spin 
effects in Mott insulators, and to the search for the mecha- 
nism of high-Tc superconductivity. To lowest order, the un- 
doped cuprate superconductors can be described by the canon- 
ical spin 1/2 two-dimensional (2D) square lattice nearest- 
neighbor (NN) Heisenberg antiferromagnet, which is among 
the simplest and most studied models in magnetism [1[]. The 
ground state displays weak classical order, reduced by quan- 
tum fluctuations at zero temperature and destroyed by ther- 
mal fluctuations at finite temperature. A possible cross-over 
between renormalized classical |2] and quantum critical |3] 
scaling, was tested experimentally in the undoped cuprates 
Sr2Cu02Cl2 (SCOC) K] and La2Cu04 (LCO) |5], and in the 
organometaflic salt Cu(DCOO)2-4D20 (CFDT) |6]. How- 
ever, while the latter shows only nearest neighbor coupling, 
high-energy inelastic neutron scattering (INS) data on LCO 
13] suggest that further-neighbor magnetic interactions influ- 
ence the above mentioned scaling measurements. Fustrated 
further-neighbor interactions could also bring the undoped 
cuprates closer to the valence bond liquid proposed as mech- 
anism for superconductivity [8]. 

It is therefore timely to investigate the excitation spectrum 
of SCOC, as an important model system. Inelastic neutron 
scattering (INS) measurements of SCOC have been limited to 
low energies and small momenta around the ordering wave 
vector 1141]. In this Letter we report the full magnetic exci- 
tation spectrum measured by resonant inelastic x-ray scatter- 
ing (RIXS). We discover a surprisingly large dispersion along 
the magnetic Brillouin zone (MBZ) boundary. An analysis 
of the data in terms of an extended Hubbard model yields a 
quantitative estimation of sizeable further-neighbor electronic 
hopping. The resulting series of longer-ranged magnetic inter- 
actions enhance quantum fluctuations, in agreement with the 
reduced ordered moment. The importance of quantum fluc- 
tuations is further revealed by differences in the spectral line 
shapes at the (-7r,0) and (-7r/2,7r/2) MZB points. 

SCOC is an insulating single-layer parent compound of the 
high-Tc superconducting (SC) materials. It is isostructural to 



the high-temperature tetragonal phase of LCO, with La re- 
placed by Sr, and apical oxygens replaced by CI. The distance 
between adjacent Cu02 planes is 18% larger than in LCO. 
AFM order develops below TAr=256 K with reduced (0.34 
IIb) moments aligned along the (110) direction in the Cu02 
plane. Both the in-plane (XY) anisotropy and the interlayer 
coupling are very small, and SCOC is an almost ideal realiza- 
tion of an 5* = 1/2 2D square-lattice Heisenberg AFM 
Magnetic excitations in SCOC have been identified by optical 
spectroscopies. A structure at 0.35 eV in absorption l^J is in- 
terpreted as a quasi-bound state of two magnons assisted by 
a phonon with momentum Qph ^ (7i',0) ]10]. Two-magnon 
excitations at E2a/=0.35 eV in the Raman spectra suggest a 
superexchange energy J ^ 0.13 eV iflTi [T2I1 . Neither Ra- 
man nor optics are sensitive to single magnons, which are 
optically-forbidden spin-flip (AS" = 1) excitations. 

RIXS in transition metal (TM) oxides probes excitations 
with mixed c harg e and spin character across or inside the 
Mott-gap ifTsUBl] . At the TM L2,3 {2p 3d) edges the 
large 2p spin-orbit interaction couples the angular momen- 
tum of the photon to the electron spin, and pure spin-flip 
excitations are possible lfT6l [TtIi . In the cuprates, dispersive 
AS = O excitations have been observed in LCO both at the 
Cu K (Is) 1 1 8,] and L3 1, 19] edges, and in the ladder compound 
Sri4Cu2404i ||20|] . More recent theoretical and experimental 
work has shown that L2,3 edge RIXS can be used to map the 
dispersion of single magnons in 2D cuprates ]21]. Therefore 
x-rays are an interesting alternative to neutrons, with the ad- 
vantage that RIXS can be performed also on sub-mm'^ sam- 
ples. 

Measurements were performed at the SAXES end station 
of the ADRESS beam line of the Swiss Light Source (SLS) 
||22|] . Single crystals (4x4x0.5 mm'^) grown from the flux 
with the c-axis perpendicular to the large surface, and charac- 
terized by x-ray and neutron diffraction, were mounted on a 
flow cryostat, with the c-axis and either the (100) or the (110) 
directions in the horizontal scattering plane. By adjusting the 
undulator, data were taken with incoming polarization either 
perpendicular (a) or within (tt) this plane. At the fixed scatter- 
ing angle of 130° the transferred momentum was ||Q||=0.85 
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FIG. 1 : (color online) (a) Schematics of the scattering geometry, (b) 
Cu L3 RIXS spectra (T=15 K) along the (100) direction. The inci- 
dent energy is set at the maximum of the absorption (XAS) (inset), 
(c) Intensity map extracted from (b). The red line is the spin-wave 
dispersion for the NN Heisenberg model and J=130 meV. 

A^^. By a rotation around a vertical axis, the projection Q|| 
of Q on the ab plane was varied in the range ± 0.73 A^^ (± 
0.92 7r/a). Positive (negative) Q values correspond to a grazing 
emission (incidence) geometry. The combined energy resolu- 
tion was AE=130 meV, and the accuracy on the energy zero 
was ± 7 meV, as determined from the elastic peak measured 
on a co-planar polycrystalline carbon sample. The momentum 
resolution, determined by the detector size, was better than 
±2.5xlO^'^(7r/a). Post-cleaving in situ or in air produced 
very similar results, confirming the good bulk sensitivity of 
L-edge RIXS. 

Figure lb is an overview of the spectra for the (100) di- 
rection for TT-polarization They were normalized to the same 
integrated intensity in the 1-2.5 eV energy range, to remove 
intensity variations due to the angular dependence of absorp- 
tion, and any other angular or time dependence. We verified 
that self-absorption, i.e. the partial re-absorption of the scat- 
tered beam, does not affect the determination of excitation en- 
ergies over the whole range of Fig. 1. The main feature at 
~1.5 eV is the manifold of optically-forbidden dd electron- 
hole excitations lfT5ll23ll . In the 0-1.5 eV energy range, where 
no electronic excitations are expected, the spectra exhibit a 
loss feature dispersing symmetrically from Q||=0 (Fig. Ic). 
Near the zone boundary, spectral weight extends beyond 300 
meV, well above the highest phonon mode (70 me V) in SCOC 
ll24ll . Its maximum follows the calculated spin- wave disper- 
sion for J=130 meV (red line, see below), strongly suggesting 
a magnetic origin of this feature. 

A model independent analysis was performed by fitting the 
main peak of the spectrum to a resolution-limited gaussian 
line shape representing the single-magnon (M) contribution 
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FIG. 2: (color online) (a) (top) RIXS spectra for a and tt polarization 
and Q||=0.58 along (100). The dashed line is the difference be- 
tween the TT spectrum and the gaussian line shape (M, solid blue line) 
representing the single-magnon contribution, (bottom) Comparison 
of Cu L-edge and O K-edge RIXS for Qy =0. 19 A'V (b) RIXS spec- 
tra at 0.92(— TT, 0) and (— 7r/2, n/2) measured with a polarization. 
The two magnons (2M) to single magnon (M) intensity ratio is 0.26 
at (-7r/2, 7r/2), and increases to 0.49 at 0.92(-7r, 0). 



(Fig. 2 top). Subtracting this line shape from the raw spectrum 
yields asymmetric features (dashed line) on both sides of the 
magnon peak. The low-energy one at ^ 50 meV, contains the 
elastic line, phonon losses and possibly a residue due to the 
approximate magnon line shape. The higher-energy feature 
reflects two-magnon (2M) and higher-order excitations, which 
give rise to the continuum above the single-magnon dispersion 
curve in Fig. Ic. Details of this analysis have negligible im- 
pact on the extracted magnon energy. The typical uncertainty 
is ±10 meV. For LCO, a similar approach yields a magnon 
dispersion in excellent agreement with INS |21]. 

An interesting and unexplored aspect of RIXS is the pos- 
sibility to separate single- and multi-magnon contributions 
by exploiting their different edge- and polarization-dependent 
cross sections, as illustrated in Fig. 2 (a). The top panel shows 
spectra (Qy — 0.58 A^^) for two polarizations, normalized 
in the 2M region for ease of comparison. The relative inten- 
sity of the single-magnon peak is reduced for a polarization. 
The bottom panel compares Cu L3 and O K edge (Is, 530 eV; 
AE=60 meV) data, for the same Q|| (0.19 A'^). Similar to 
the Cu K-edge case, the O K edge spectrum contains only the 
2M continuum, which was suggested by previous experiments 
lHHH^, and is now clearly resolved around 0.4 eV for this 
Q||. By contrast, the L3 line shape exhibits a prominent sin- 
gle magnon loss at ^ 0. 1 e V. Figure 2 (b) ilustrates the spec- 
tral line shapes near the zone boundary points (— 7r/2, 7r/2) 
and {—TT, 0). Near (— tt, 0) the magnon peak is weaker, and 
spectral weight is transferred to the higher energy continuum. 
The comparison strongly suggests that the quantum effect ob- 
served by neutrons in CFTD 12711 is also present in SCOC, and 
most likely in the cuprates in general. This observation is im- 
portant for theories arguing that superconductivity occur from 
strong magnetic fluctuations |28, 29], although further devel- 
opment of the RIXS cross-section and more quantitative pre- 
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FIG. 3: (color online) a) Magnon energies extracted from the RIXS 
data. Open and closed symbols stem from 2 independent measure- 
ments on different samples. Dot-dashed line is a NN Heisenberg 
model with J = 130 meV. Red dashed line is a NN Hubbard model 
fit for t = 0.261 ± 0.004 eV and U = 1.59 ± 0.04 eV. Blue lines 
are the further neighbor Hubbard fits, b) same as a) with NN Heisen- 
berg dispersion subtracted to better visualize details of the dispersion. 
The blue band shows the spread in dispersions obtained for fits with 
1.9 eV< t/ < 4 eV. c) The fit to the neutron data on La2Cu04 |7], 
shifted by (vr, tt) for ease of comparison. 
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FIG. 4: (color online) Results of extended Hubbard model fits for 
SCOC (left) and LCO (right), a-b) after fitting t and t' for fixed 
U and t" . Dashed lines mark the minima versus U , giving c-d) and 
e-f). c-d) versus IJ , showing the range of U that our data can 
support, e-f) extracted hopping parameters t, t' and t" as function 
of U. For given U, the uncertainty on the ts is around 5 meV. In all 
panels, thick red lines are for t" > and thin black lines for t" < 0. 



dictions from these theories are required before definite con- 
clusions can be drawn. 

The Q-dependence of the magnon energy extracted from 
the data of Fig. 1, and from similar data for the (110) di- 
rection and two different samples, is summarized in Fig. 3. 
The RIXS data are consistent with the small-Q results from 
INS, but cover for the first time the full dispersion up to the 
boundary of the MBZ. They reveal a striking 70 meV differ- 
ence between the magnon energies of 310 meV at (7r,0) and 
240 meV at (7r/2,7r/2). This can be compared with the smaller 
^-^20 meV dispersion in LCO [7]. Dispersion along the zone 
boundary in all cuprates is also predicted by recent theory, 
which however underestimates the (7r,0) energy in SCOC by 
almost 50 meV [M- 

For the simple 5*= 1/2 2D Heisenberg model with NN ex- 
change, linear spin-wave theory predicts a constant magnon 
energy fujj=2J along the MBZ boundary. First order quan- 
tum corrections uniformly renormalize the dispersion by a 
factor Zc=1.18. Numerical results llsil l32ll and neutron data 
on CFDT iIztI [33ll have established that the magnon energy 
for purely NN exchange is actually 6% larger at (7r/2,7r/2) 
than at (7r,0). The dispersion in Fig. 3 and in LCO is in the 
opposite direction to this quantum effect. It could be repro- 
duced by adding freely adjustable further-neighbor exchange 
interactions. However, the Heisenberg (spin-only) hamilto- 
nian is the low-energy projection of an electronic system at 
half filling, and a better approach is to systematically con- 
sider higher orders to this projection. Indeed the dispersion 
in LCO ^ was described by projecting the one-band Hub- 
bard model with effective Coulomb repulsion U to 4th order 
in the NN hopping t, giving rise to further-neighbor exchange 
interactions J, J' = J" and Jc- The same approach gives 
for SCOC an unphysically low value of f/ = 1.59 ± 0.04 eV 



and t = 0.261 ± 0.004 eV. This indicates that the Hubbard 
model with only NN hopping has reached its limit. A more 
plausible approach is to include further-neighbor hoppings t' 
and t" 1 34]. We therefore extended the analysis to 4th order 
in t, t' and t" and find that this approach give a more reason- 
able range of U , is consistent with higher energy RIXS and 
ARPES, and provide better fits to the shape of the dispersion 
(Fig. 3b). 

Magnetic excitations provide accurate information on the 
interactions, but do not directly probe U . A more direct 
probe of the effective U was found in the higher energy 
part of the RIXS Cu-K spectra from the sister compound 
Ca2Cu02Cl2, where a 2.5-4 eV dispersive feature analyzed 
within the one-band Hubbard model was reported consistent 
with U = 3.5 eV 1.141 . However, both the experimental and 
the numerical accuracy of such determinations could be im- 
proved in the future. Therefore we performed fits as a function 
of fixed U, as summarized in Fig. |4] The dispersion is sym- 
metric in the signs of t' and t", and leads to two possible solu- 
tions with t't" < and t't" > respectively. In compliance 
with ARPES and theoretical estimates, we assume t' < 0, and 
lean towards t" > 0. Both solutions for t"t ^ constrain U 
to larger than 2 eV and smaller than 4 eV for t"t > 0, 
and give essentially the same t/U and t' /t, which depend only 
weakly upon the chosen U, respectively from 0. 17 to 0. 12 and 
from -0.31 to -0.42. Hence, our data provide strict constraints 
on the effective parameters that can be used in the one-band 
Hubbard model for SCOC: U larger than 1.9 eV; significant 
second neighbor hopping \t' /t\ > 0.31; and a unique set of 
hopping parameters for a given U (Fig. 15). For U = 3.5 eV, 
we obtain t/C/ = 0.139 ± 0.004, i'/t = -0.41 ±0.01 or -0.38 
±0.01 and t"/t = 0.14±0.01 or -0.32±0.01. These param- 
eters are roughly consistent with ARPES results from SCOC 
ll35ll . which were described hy U — 3.5 eV, t = 0.35 eV, 



4 



t' /t = —0.35 and t" /t = 0.22, with the accuracy of the com- 
parison hkely set by the broad ARPES Hnewidth. Thereby, we 
have derived a consistent description within a single model of 
both spin-wave and ARPES spectra. 

To gain insight into the origin of the zone-boundary disper- 
sion, we performed the same extended Hubbard model anal- 
ysis for LCO (Fig.lHright), using the 10 K data by Coldea et 
al. 10] • Also for LCO we find that the spin-wave dispersion 
can be described by a larger U more compatible with higher 
energy probes, and that there is significant further neighbor 
hopping, albeit slightly weaker than in SCOC. This trend, and 
our quantitative hopping parameters, are consistent with re- 
cent LDA and LDAh-U calculations ll36l [37ll . confirming ex- 
perimentally that the main effect of the apical atom (oxygen in 
LSCO, chlorine in SCOC) is to influence the further-neighbor 
hopping terms t' and t" . For the projected Heisenberg hamil- 
tonian, this means that the difference in spin-wave dispersion 
comes not from different main interactions J and Jc, but from 
the many additional further-neighbor hopping paths and hence 
magnetic interactions. The existence of so-called 'ring ex- 
change' Jc coupling 4 spins on a plaquette generated signifi- 
cant theoretical efforts to understand the consequences of this 
term. An important conclusion from our work is that this term 
is not unique and that many further interactions and larger 4- 



spin loops have similar weight. 

In summary, we have employed RIXS to obtain qualitative 
and quantitative new insight into the magnetic excitation spec- 
trum of the representative two-dimensional antiferromagnetic 
insulator Sr2Cu02Cl2. Measuring the full spin-wave spec- 
trum, we found a large zone boundary dispersion, which we 
could reproduce with an extended t-t'-t"-U Hubbard model, 
placing quantitative constraints on the hopping parameters. 
Most notably, we demonstrate sizeable longer range hopping, 
and henceforth magnetic interactions. Taking into account the 
quantum corrections generated by these higher-order hopping 
terms is essential to achieve a quantitative description of the 
ground state properties of SCOC, namely the reduced value 
of its ordered moment lf34ll . In a broader perspective, these 
results establish an important reference and suggest a general 
method, requiring only small crystals, to address the nature of 
the ground state, and the evolution of magnetic correlations 
throughout the phase diagram of the cuprates. 

The SAXES insti'ument at the ADRESS beamline of the 
Swiss Light Source was jointly built by Paul Scherrer Institut 
and Politecnico di Milano. We gratefully acknowledge dis- 
cussions with M. Mourigal, M. Gingras, J.-Y. Delannoy, F. 
Vernay and B. Normand, and support from the Swiss NSF 
and the MaNEP NCCR. 
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